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abstract 

Thirty-four  near-simultaneous  pairs  of  CTD  and  Sippican  mrxlel  T-5  XBT  profiles  were  obtained  during  an 
experiment  in  the  Sargasso  Sea  during  the  summer  of  1991.  The  data  were  analyzed  to  assess  the  temperature 
and  fall-rate  accuracies  of  the  T-5  probes.  The  XBT  temperatures  averaged  O  OVT  warmer  than  CTD  temper¬ 
atures,  with  some  suggestion  that  the  offset  might  be  different  for  different  acquisition  systems  and  that  it  might 
be  slightly  temperature  or  pressure  dependent.  When  the  offset  was  removed,  the  differences  betw<*»n  CTD  and 
XBT  temperatures  had  a  standard  deviation  of  about  0.08°C  over  a  temperature  range  of  .T’-29°f' 

An  improved  elapsed  fall-time-to-depth  conversion  equation  for  Sippican  I-5's  in  the  Saigasso  Sea  was  found 
to  be  c  =  6.705/  -  0.001619/’.  with  r  the  depth  in  meters  and  /  the  elapsed  fall  time  of  the  probe  in  seconds 
The  standard  deviation  of  depth  was  about  8  m  over  a  depth  range  of  0  to  approximately  1800  m.  A  cubic  fit 
to  the  data  was  equally  good  or  slightly  better.  Whether  a  geographically  universal  fall-rate  equation  can  be 
devised  for  each  model  XBT  is  still  unclear.  In  addition,  now  that  a  number  of  different  manufacturers  are 
introducing  their  versions  of  XBTs  and  XBT  acquisition  systems  onto  the  market,  unresolved  questions  exist 
regarding  the  differences  in  data  taken  with  these  different  models. 


1.  Introduction 

The  ship-deployed  expendable  bathythermograph 
(XBT)  is  a  relatively  inexpensive  nonrecoverable  in¬ 
strument  that  yields  temperature-versus-depth  values 
down  to  a  maximum  of  about  1 800  m,  depending  upon 
the  model.  XBTs  were  first  introduced  in  1965  by  Sip¬ 
pican  Ocean  Systems,  and  over  four  million  have  been 
sold.  Recently,  other  manufacturers  have  also  begun 
to  enter  the  market.  Global  archives  are  dominated  by 
XBT  measurements,  and  many  oceanographic  exper¬ 
iments  depend  upon  them  as  the  primary  tool  for  ob¬ 
taining  subsurface  measurements.  The  XBT  is  truly 
the  “work  horse”  of  physical  oceanography.  Clearly, 
knowledge  of  the  accuracy  and  possible  biases  of  this 
instrument  is  very  important. 

A  quick  description  of  the  generic  XBT  will  help 
clarify  where  some  of  the  inaccuracies  and  errors  in 
measurements  might  arise.  A  thermistor  is  mounted 
inside  the  nose  of  a  small,  streamlined,  weighted  body, 
which  upon  deployment  is  released  from  its  surround¬ 
ing  canister  and  falls  over  the  side  of  a  ship.  As  the 
probe  descends,  a  thin,  two-stranded  insulated  wire 
unwinds  from  two  spools,  one  located  in  the  falling 
body  and  the  other  in  the  canister  that  remains  topside 
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on  board  the  ship.  The  thermistor's  resistance  changes 
as  the  probe  falls  through  water  of  different  tempera¬ 
tures,  and  the  varying  voltage  required  to  maintain  a 
constant  amperage  is  sensed  topside.  This  variation 
may  be  used  to  generate  a  trace  on  a  special  chan  re¬ 
corder  or  may  be  translated  by  an  analog-lo-digital  (A/ 
D)  converter  into  a  digital  reading.  These  readings  are 
then  translated  into  temperature  values.  Depth  is  cal¬ 
culated  from  the  assumed  fall  rate  of  the  weighted  body  , 
it  is  not  measured  directly.  For  its  XBTs,  Sippican 
Ocean  Systems  claims  an  overall  temperature  accuracy 
of  0. 15°C  and  a  depth  accuracy  of  ±5  m  or  2T  of 
depth,  whichever  is  larger  (Sippican  Ocean  Systems 
1991;  Sippican  Ocean  Systems  1992,  personal  com¬ 
munication). 

Temperature  inaccuracies  and  biases  may  originate 
I )  in  the  thermistor  itself,  2 )  in  transient  responses  due 
to  the  thermistor  time  constant  (and  possibly  other 
electronic  components)  and  to  the  thermal  inertia  of 
the  probe,  3)  in  the  voltage  changes  that  may  be  con¬ 
taminated  by  mechanical  and  electrical  interference 
and  by  damage  to  the  integrity  of  the  wire  insulation, 
and  4)  through  errors  in  converting  the  topside  voltage 
into  either  a  strip-chart  temperature  or  a  digitally  re¬ 
corded  temperature.  A  detailed  analysis  of  many  of 
these  problems  and  how  they  are  manifested  in  the 
temperature  trace  have  been  described  by  Blumenthal 
and  Kroner  (  1978).  Roemmich  and  Cornuelle  (1987) 
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examined  the  ultimate  accuracy  that  might  be  achievc-d 
with  conventional  XBTs  with  frequent  calibration  of 
an  electronic  digitization  system  and  predeployment 
calibration  of  each  probe  itself.  They  coni'luded  that 
total  temperature  error  could  be  reduced  to().()2°('  or 
less,  but  depth  errors  due  to  fall-rate  inaccuracy  could 
degrade  this  to  0. 1  °-().5°C. 

Depth  accuracy  or  lack  thereof  is  due  directly  to  the 
suitability  of  the  fall-rate  equation  and  to  the  effect  of 
transients  caused  by  the  initial  conditions  of  the  probe 
when  it  enters  the  water.  Theoretical  approaches  toward 
understanding  the  .XBT's  fall  rate  have  been  given  by 
Green  ( 1984)  and  Hallock  and  Teague  ( 1992). 

A  number  of  observational  studies  of  the  achieved 
temperature  and  depth  accuracies  of  various  types  of 
.XBTs  have  been  published.  Most  have  e.xamined 
probes  made  by  Sippican  Ocean  Systems  [Hanawa  and 
Yoritaka  (  1987)  being  an  exception],  and  most  have 
been  T-7's,  which  have  a  maximum  depth  range  of 
about  760  m.  One  of  the  first  published  papers  w-as  by 
Rierl  and  Robinson  (  1977).  who  noted  both  random 
and  systematic  differences  in  isotherm  depths  derived 
from  simultaneous  CTD  (conductivity,  temperature, 
and  depth)  and  Sippican  T-7  XBT  data  in  the  North 
Atlantic.  Random  errors  were  on  the  order  of  8  db 
(about  8  m).  while  systematic  differences  suggested  the 
XBTs  fell  more  rapidly  than  expected  down  to  about 
400  m.  and  less  rapidly  after  that.  Heinmiller  et  al. 
(198.2)  examined  Sippican  T-4  (~450  m)  and  T-7 
profiles  as  compared  with  simultaneous  CTD  casts  in 
the  tropical  Pacific.  They  found  the  XBTs  to  be  con¬ 
sistently  warmer  than  the  CTDs,  and  after  correcting 
the  temperatures,  they  found  depth-dependent  differ¬ 
ences  between  CTD  and  XBT  isotherm  depths,  with 
the  nominal  XBT  depth  being  consistently  too  shallow 
below  several  hundred  meters.  Both  of  these  studies 
used  XBT  paper  traces,  which  may  be  characterized 
by  somewhat  different  errors  than  digitally  recorded 
data. 

Hanawa  and  Yoritaka  {  1987)  compared  CTD  data 
from  the  northwest  Pacific  with  that  from  T-7  equiv¬ 
alent  XBTs  made  by  a  Japanese  manufacturer  and  us¬ 
ing  a  digital  XBT  system.  They  found  lot-dependent 
temperature  biases  and  proposed  a  corrected-depth 
equation.  This  work  was  followed  by  Hanawa  and 
Yoshikawa  ( 1991 ),  who  examined  the  results  of  five 
CTD-Japanese  T-7  XBT  comparison  experiments,  all 
in  the  northwest  Pacific  but  in  locales  with  quite  dif¬ 
ferent  temperature  gradients.  They  concluded  that  the 
nominal  fall-rate  equation  yielded  an  underestimate  of 
the  true  depth  over  much  of  the  profile,  and  devised  a 
revised  equation.  They  found  that  different  datasets 
yielded  different  empirical  equations;  however,  they 
did  not  present  error  bars  on  the  empirical  equation 
coefficients  to  show  that  the  equations  were  truly  sta¬ 
tistically  different.  Their  findings  may  be  relevant  to 
T-7  probes  made  by  other  manufacturers,  but  this 
needs  verification. 


Wright  and  S/abados  (  1989  )  examined  the  temper¬ 
ature  and  depth  accuracies  of  several  types  of  Sippican 
.XBTs  in  a  part  of  the  northw  est  tropical  Atlantic  w  here 
salt  fingering  leads  to  thermohaline  staircases,  which 
are  sections  of  extremely  well-mixed  layers  ranging 
from  several  meters  to  several  tens  of  meters  in  thiek- 
ne.ss  separated  by  thin  ( several  meters  thick )  interfaces. 
Their  data  came  from  modest  numbers  of  1-4.  15. 
T-6,  T-7.  and  T-10  probes  and  four  different  acquisition 
systems.  The  extremely  well-mixed  layers  allowed  them 
to  make  a  very  good  estimate  of  the  achievable  tem¬ 
perature  accuracy  of  the  probes  and  the  determination 
of  any  sort  of  bias.  All  versions  of  the  probes  seemed 
to  give  temperature  values  a  bit  warmer  than  the  C  TD 
values,  and  these  biases  varied  from  system  to  system, 
but  due  to  the  small  sample  sizes  (5-15  probes),  the 
statistical  significance  of  these  differences  is  equivocal. 
Particularly  relevant  to  this  study,  they  found  a  T-5 
temperature  bias  of  0.11  °C  (standard  deviation  of 
0.06°C)  for  a  Sippican  MK-9  acqufsition  system  and 
0.24°C  (0. 17°C)  for  a  Bathysystems  unit.  For  all  but 
the  T-5  probes  their  findings  were  similar  to  those  of 
other  researchers  in  indicating  that  the  actual  fall  rates 
were  faster  than  the  nominal  fall  rates.  The  T-5.  how¬ 
ever,  seemed  to  fall  somewhat  slower  than  the  nominal 
rate. 

Singer  ( 1990)  compared  CTD  and  Sippican  T-7  data 
from  the  Gulf  of  Mexico  by  comparing  depths  of  iso¬ 
therms.  Two  different  manufacturers'  acquisition  sys¬ 
tems  were  used,  and  no  attempt  was  made  to  check 
for  possible  temperature  biases.  His  findings  were 
broadly  in  agreement  with  other  work  in  that  he  found 
that  below  about  lOO  m  the  nominal  fall-rate  equation 
yielded  depths  that  were  shallow  compared  with  the 
CTD  depths. 

As  this  review  indicates,  most  past  studies  have  con¬ 
centrated  on  the  characteristics  of  T-7  XBTs.  The  T-7 
temperatures  appear  to  be  slightly  higher  than  CTD 
temperatures,  and  the  actual  fall  rates  appear  to  be 
somewhat  faster  than  the  nominal  fall  rate  supplied  by 
the  manufacturers.  In  our  study,  we  examined  34  T-5 
profiles  and  compared  them  with  25  concurrent  CTD 
profiles.  In  section  2,  we  describe  the  XBT  and  CTD 
datasets  and  their  processing.  In  .section  3.  wc  pre.sent 
our  analyses  of  the  temperature  and  depth  errors  in 
the  XBT  datasets,  and  in  section  4,  we  summarize  our 
findings  and  discuss  their  implications. 

2.  Data  sources  and  processing 

During  an  experiment  in  the  Sargas.so  .Sea  south  of 
Bermuda  during  summer  1991. 34  Sippican  T-5  XBTs 
were  dropped  nearly  simultaneously  with  25  CTD  casts 
in  order  to  develop  improved  T-5  fall-rate  equations 
for  the  area  and  to  assess  XBT  temperature  errors  or 
biases.  Three  different  ships  were  involved,  with  three 
different  CTD  systems  and  three  different  XBT  sys¬ 
tems.  as  summarized  in  Tabic  1 .  The  CTD  casts  were 
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lengthy,  going  to  at  least  2000  m,  and  many  to  over 
5000  in.  Hence,  “simultaneous"  ended  up  meaning 
within  8  km  and  3.25  h  of  the  nominal  start  of  the 
CTD  cast  ( Fig.  I ),  Larger  space  and  time  windows  led 
to  anomalously  large  deviations  between  features  mea¬ 
sured  by  the  CTDs  and  corresponding  features  mea¬ 
sured  by  the  XBTs.  As  it  is.  there  is  a  certain  amount 
of  inherent  uncertainty  in  the  results  that  is  due  to 
oceanic,  not  XBT,  variability.  This  is  probably  on  the 
order  of  ±5  cm  in  depth  (mostly  internal  wave  in¬ 
duced  )  and  ±0.05°  to  ±0, 1  °C  in  temperature  ( assum¬ 
ing  a  temperature  gradient  of  0.01°-0.02°C  m  '). 


Time  (Mrs) 

Flo.  1.  Temporal  and  spatial  separations  between  the  nominal 
start  of  the  CTD  ca.sts  and  the  deployment  of  their  a'  .oeiated  XBTs, 


Since  this  variability  is  probably  randomly  distributed, 
it  should  not  add  a  bias  to  the  results,  but  it  will  increase 
the  achievable  error  bounds  on  our  estimates. 

Standard  processing  procedures  produced  temper¬ 
ature  ( and  salinity )  versus  pressure  profiles  for  the  C  TD 
data  from  all  three  ships.  All  C  LDs  had  been  recently 
calibrated,  although  the  second  Lynch  cast  was  dtine 
using  the  backup  CTD  w  hose  last  calibration  dale  was 
uncertain  but  certainly  within  the  preuous  year.  After 
considering  the  results  of  the  recent  calibrations,  we 
consider  the  CTD  temperatures  to  be  accurate  to  at 
least  ±0.()()5°C.  and  pressures  to  ±4.5  db.  During  the 
proce.ssing  the  samples  were  "bin  averaged"  into  1-db 
bins.  Pressure  was  then  converted  to  depth  using  the 
methtxl  described  in  Saunders  (  1981 ).  anu  the  resulting 
profiles  were  intcrjyolated  { linearly )  to  a  i-m  depth  in¬ 
terval. 

Nominal  XBT  temperatures  r  .d  depths  were  cal¬ 
culated  by  either  the  Sippican-s  .pplied  software  in  the 
case  of  the  acquisition  sy, stems  on  board  the  T  7  and 
Ran^c  Rover  or  the  NA  VOCEANO-supplied  software 
for  the  system  on  board  the  Lynch  { Table  I ).  .Ml  sys¬ 
tems  calculated  nomir.  l  XBT  depths  from  the  elapsed 
fall  time  of  the  prolx  using  the  manufacturer's  supplied 
equation 

-  6.828/  -  0.()()I82/n 

where  r  is  d.;pth  in  meters  and  /  is  elapsed  fall  time  in 
seconds,  '-.lapsed  fall  time  was  later  backed  out  from 
the  depth  by  solving  the  above  quadratic  for  /  and 
picking  the  appropriate  root,  which  is  unambiguous. 
The  XBT  sampling  rate  was  10  Hz.  leading  to  a  spacing 
of  approximately  64  cm.  so  the  profiles  were  also  lin¬ 
early  interpolated  to  a  I-m  spacing. 

Corresponding  features  on  the  simultaneous  XBT 
and  CTD  profiles  were  matched.  Since  much  of  the 
smalier-scale  structure  in  an  ocean  profile  is  masked 
by  the  large-scale  structure  of  the  thcrmocline,  we 
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FiCi.  2.  Example  of  an  unfiltered  CTD  temperature  profile  from  this  studs,  its  bandpass-filtered 
version,  and  the  bandpass-filtered  version  of  itsas.sociated  .XBT  drop.  Matched  features  arc  indicated 
on  the  filtered  profiles 


bandpass-filtered  both  CTD  and  XBT  profiles  with  a 
simple  boxcar  filter,  a  technique  suggested  to  us  by  M. 
Prater  and  later  described  in  Prater  (1991).  Half-power 
points  were  chosen  at  5  and  100  m.  An  example  of  this 
matching  of  filtered  profiles  is  illustrated  in  Fig.  2. 
Features  were  chosen  between  approximately  20  and 
1 775  m,  and.  as  Fig.  2  indicates,  they  were  distributed 
as  evenly  over  the  full  water  column  as  possible.  Ap¬ 
proximately  16  points  per  profile  were  selected.  The 
CTD  and  nominal  XBT  depths  at  which  the  features 
matched  were  recorded,  along  with  the  imfiltered  tem¬ 
perature  values  at  those  depths.  The  result  was  a  set  of 
5-14  observations,  each  observation  consisting  of  CTD 
depth.  CTD  temperature.  XBT  nominal  depth,  XBT 
elapsed  fall  time,  and  XBT  temperature.  Compari.son 
of  the  boxcar  filter  results  with  other  possible  filters 
indicated  little  if  any  difference  in  the  resulting  set  of 
observations.  Because  corre.sponding  features  were 
matched,  not  isotherms,  the  fall-rate  equation  findings 
are  not  contaminated  by  errors  due  to  temperature  er¬ 
rors,  as  would  be  the  case  if  isotherms  were  matched. 

Prater  ( 1991 )  goes  on  to  describe  an  intuitively  ap¬ 
pealing  correlation  procedure  for  shifting  vertical  sec¬ 
tions  of  corresponding  filtered  profiles  up  and  down  in 
depth  until  the  correlation  is  maximized  and  for 
choosing  the  depth  offset  at  which  the  correlation  is 
maximized  as  the  depth  ofl'set  between  the  CTD  and 
XBT  depths.  We  tried  this  procedure  but  found  it  gave 
a  considerably  larger  ambiguity  in  depth  than  matching 
the  profiles  manually.  Perhaps  the  algorithm  can  be 
refined  to  eliminate  this  difference,  but  wc  opted  to  use 
the  straightforward  direct  measurement  approach. 


3.  Results 

a  Tenipcraluri’  accuracy 

The  CTD-XBT  feature  temperature  differences  arc 
plotted  versus  temperature  in  Fig.  3a.  .A  small  tem¬ 
perature  correction  appears  to  be  needed,  which  might 
be  a  function  of  true  temperature  or  of  depth  ( probably 
pressure,  in  actuality).  Correcting  by  temperature  is 
more  useful,  since  it  y  ields  equations  of  the  form  7',.an 
=  «  -f  hJ'xBT  relating  corrected  XBT  temperature  to 
measured  XBT  temperature  (Table  2).  One  of  the  ac¬ 
quisition  systems  {E  E  )  had  a  statistically  significant 
dependence  on  /"xb!  '.  the  other  two  ( l.yncli  and  Ranyc 
Rover)  did  not.  and  for  these  data,  simply  correcting 
with  an  offset  would  be  justifiable.  For  consistency  we 
used  the  linear  equations  for  all  datasets.  The  corrected 
temperatures  are  shown  (in  Fig.  lb)  to  be  significant 
improvements  over  the  uncorrected  temperatures. 

Kennelly  et  al.  ( 1989)  and  Wright  and  Szabados 
( 1989)  also  reported  temperature  biases  with  Sippican 
T-5  XBTs.  Kennellv  el  al.  (  1989)  reported  a  bias  of 
0.075“C.  and  Wright  and  Szabados  (  1989)  of  0.1 1° 
and  ().24'’C  (different  values  for  different  acquisition 
systems),  with  the  XBT  temperatures  being  warmer 
than  the  CTD  temperatures.  T  hese  arc  of  the  same 
magnitude  as  our  mean  CTD- .XBT  differences  ( Table 
2).  although  perhaps  somewhat  larger.  While  we  did 
not  record  the  lot  numbers  of  our  XBTs.  ours  almost 
certainly  did  not  come  from  the  same  lots  as  the  other 
two  studies,  since  our  probes  had  been  obtained  at  least 
two  to  three  years  later  than  the  probes  in  the  earlier 
studies.  After  correcting  the  temperatures,  the  pooled 
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Fk).  3,  (a)  The  difference  berween  CTD  temperature  minus  XBT  uncorrected  temperature  plotted  versus  ('TD  depth  for  all  CTD-XBl 
pairs.  The  positive  and  negative  S'*  of  the  outliers  have  been  trimmed.  Here  ■  indicates  Lynch.  A  indicates  I-  I'  and  •  indicates  Run^c 
Rover  (b)  Same  as  (a)  except  for  corrected  temperatures. 


Standard  deviation  of  our  XBT  temperatures — repre¬ 
senting  thermistor  to  thermistor  variability  and  the  ef¬ 
fects  of  small-scale  changes  in  the  ocean  thermal  struc¬ 
ture  between  CTD  and  XBT  measurements — was 
0.08°C.  Kennelly  et  al.  ( 1989)  reported  an  rms  vari¬ 
ation  of  less  than  0. 1  °C,  and  Wright  and  Szabados 
( 1989)  a  standard  deviation  of  0.06°  and  0.1 7°C. 

h.  Fall-rale  equation 

The  usually  assumed  form  of  the  fall-rate  equation 
is  a  quadratic  forced  through  zero,  that  is,  an  equation 

Table  2.  Statistics  of  the  differences  between  the  CTD  feature 
temperatures  and  the  XBT  feature  temperatures,  and  the  derived 
correction  equations.  The  corrected  temperature  is  of  the  form 
--  a  -E  bTxev-  The  asterisk  indicates  those  coefficients  no,  significantly 
different  from  I .  After  correction,  the  mean  differences  were  less  than 
I  X  10  ’.  Calculations  were  done  on  each  platform’s  dataset 
individually  and  then  on  the  pooled  dataset.  A  10%  trim  was  applied 
to  each  dataset  to  remove  outliers. 


Platform 

Number 

Mean  T 
difference 
before 
correction 

Correction 

coefficients 

and 

values 

Standard 

deviation 

after 

correction 

Lynch 

77 

-0.104 

a  =  -0.08771 
b  =  0.99869* 

0.069 

Range 

Rover 

127 

-0.077 

a  =  -0.06633 
b  =  0.99920* 

0.092 

ET 

277 

-0.052 

a  -  -0.02038 
h  =  0.99760 

0.070 

Pooled 

481 

-0.067 

a  -  -0.04193 
b  =  0.99811 

0.079 

of  the  form  z  =  at  -l-  hi",  with  l  the  elapsed  fall  time. 
For  the  XBTs  used  in  this  study,  the  manufacturer-— 
Sippican  Ocean  Systems — supplies  the  equation 

z  =  6.828f  -  O.OOI82f‘’. 

The  differences  between  CTD  and  XBT  feature  depths 
from  this  equation  versus  CTD  depth  are  plotted  in 
Fig,  4.  There  is  a  clear  depth-dependent  error  in  the 
nominal  XBT  depth.  To  derive  an  improved  depth 
versus  elapsed  XBT  fall-time  equation,  we  started  by 
noting  the  XBT  depth  for  each  CTD-XBT  matched 
feature  and  calculating  the  XBT  elapsed-time  value  by 
inverting  the  manufacturer’s  fall-rate  equation,  as  de¬ 
scribed  earlier.  The  CTD  depth  of  that  feature  is  con¬ 
sidered  the  “true”  depth,  with  inherent  uncertainty  due 
primarily  to  internal-wave  variability.  We  ran  several 
regressions  of  z  (true  depth)  against  t  (elapsed  fall 
time),  assuming  five  different  polynomial  regression 
models.  These  models  were  a  linear  equation,  both  with 
and  without  a  constant  term;  a  quadratic  equation, 
also  with  and  without  a  constant;  and  a  cubic.  Quartic 
and  higher  models  yielded  coefficients  for  the  highest- 
order  terms  that  were  not  significantly  different  from 
zero. 

Criteria  for  selecting  the  best  regression  model  from 
among  a  set  of  candidates  are  described  in  Kleinbaum 
et  al.  ( 1 988 ).  The  first  criterion  is  based  upon  selection 
of  the  largest  sample  squared  multiple  correlation  coef¬ 
ficient  R^.  but  it  gives  inconclusive  results,  since  all 
models  have  highly  significant  regressions  with  val¬ 
ues  of  greater  than  0.999.  The  fourth  criterion,  Mal¬ 
low’s  Cf,,  is  similarly  ambiguous,  suggesting  only  that 
more  than  a  linear  model  and  less  than  a  sixth-order 
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Fig.  4.  The  diBerence  between  CTD  depth  minus  XBT  depth 
plotted  versus  CTD  depth  for  all  CTD-XBT  pairs.  The  positive  and 
negative  S'-’r  of  the  outliers  have  been  trimmed  Here  ■  indicates 
J.ync'h.  A  indicates  £.  T .  and  •  indicates  Kanf;e  Rover. 


mcxlel  is  needed.  Criteria  2  and  3  proved  more  useful. 
Criterion  2  proposes  computing  a  test  statistic  to  com¬ 
pare  the  highest-order  model  (“maximum  model”  or 
“fe-variable  model")  with  lower-order  models.  If  the 
statistic  is  not  significant,  then  the  lower-order  model 
is  adequate.  The  test  statistic  Fp  is  calculated  according 
to 


_  [SSE(p)-SSE(k)]/(k-  p) 

”  "  MSE(/c) 

where  k  is  the  number  of  variables  in  the  highest-order 
model  (3  for  a  cubic),  p  is  the  number  of  variables  in 
the  other  models  being  considered,  SSE(/7)  is  the  error 
sum  of  squares  for  the  /j-variable  model  and  SSE(^) 
for  the  ^-variable  model,  and  MSE(^)  is  the  mean- 


square  error  for  the  A- variable  mtxlel.  This  statistic  is 
compared  to  an  /-'  distribution  with  k  -  p  and  n  -  k 
-  I  degrees  of  freedom.  The  results  are  summarized 
in  Table  3,  along  with  the  statistics  using  the  manu¬ 
facturer’s  model.  The  cubic  model  is  better  than  the 
linear  models  or  the  manufacturer's  fall-rate  mcxlel  at 
significance  levels  exceeding  99'; .  and  it  is  better  than 
the  quadratic  mtxlcls  at  a  95‘X-  level  but  not  a  99'1 
level. 

The  third  criterion  suggests  picking  the  model  that 
has  the  smallest  error  variance.  MSE.  From  Table  3 
we  see  again  that  the  cubic  model  has  a  much  lower 
MSE  than  the  linear  or  manufacturer-supplied  models 
and  a  slightly  lower  MSE  than  the  quadratic  models. 

The  quadratic  and  cubic  mcxlels  both  fit  our  data 
significantly  better  than  the  manufacturer’s  model,  with 
the  cubic  being  slightly  better.  Significance  tests  on  the 
values  of  the  coefficients  indicate  that  for  the  quadratic 
model  the  constant  coefficient  is  not  significantly  dif¬ 
ferent  from  zero.  This  leaves  the  quadratic  forced 
through  zero  and  the  cubic  models  as  the  models  from 
which  to  choose.  Both  have  residuals  almost  but  not 
perfectly  normally  distributed,  with  the  residuals  from 
the  cubic  being  slightly  closer  to  a  normal  distribution 
(Figs.  5a,b).  The  cubic  model,  then,  would  appear  to 
be  the  model  of  choice,  but  whether  or  not  it  is  suffi¬ 
ciently  better  to  justify  the  added  complexity  of  the 
cubic  term  is  moot.  Further  studies  of  this  type  with 
larger  datasets  and  improved  simultaneity  between 
CTD  and  XBT  measurements  are  needed  to  resolve 
whether  a  cubic  equation  is  more  appropriate  than  the 
usually  assumed  quadratic. 

The  results  for  the  two  best  models  are  summarized 
in  Table  4.  The  best  model  was 

c  =  -1.803  +  6.795/  -  0.002475/  = 


-I-  2.148  X  10  ^  \ 


followed  by 

z  =  6.705/  -  0.001619/=. 


Table  3.  Evaluation  of  the  candidate  fall-rate  models  according  to  criteria  2  and  3  from  Kleinbaum  et  al.  (1988).  Criterion  2  considers 
Fp,  the  test  statistic  comparing  all  other  models  with  the  cubic  model  (the  maximum  model).  Its  formula  is  given  in  the  text.  The  test  statistic 
is  compared  with  an  /•'distribution  whose  critical  values  at  a  95%  and  99%  significance  level  are  given  in  the  last  column.  Criterion  3  looks 
for  the  minimum  MSE.  The  best  models  are  the  cubic  model  and  the  quadratic  forced  through  zero. 


Model 

SSE 

MSE 

l-'r 

^  cnl 

95‘^:/99% 

z  -  hi 

0.10322  X  10“ 

0.21417  X  10’ 

487 

3.02/4.66 

e  =  a  +  br 

0.77987  X  10’ 

0.16214  X  10’ 

309 

3.02/4.66 

z  ■=  hi  +  cF 

0.34423  X  10’ 

0.71566  X  10^ 

5.3 

3.86/6.69 

z  =  a  +  ht  ^  cF 

0.34384  X  10’ 

0.71634  X  lO’ 

4.7 

3.86/6.69 

Maximum  model; 
z  =  a  +  b(  +  cri  + 

0.34047  X  10“ 

0.71080  X  10’ 

Sippican  model; 

z  =  6.828/  -0.00182r^ 

0.92001  X  10’ 

0.19127  X  10’ 

8.15  X  10* 

3.86/6.69 
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Fio.  5.  Residuals  for  the  two  best  models  found  in  this  study;  (a)  for  the  quadratic  forced  through  /cro;  ( h )  for  the  cubic  nuKiel 
Comparison  with  Fig.  A  illustrates  the  improvement  in  depth  accuracy  that  is  achieved  with  the  new  equations 


The  difference  between  the  two  equations  is  less  than 
1  m  betw'een  about  25  and  250  s.  Within  the  useful 
depth  range  of  the  T-5,  the  maximum  difference  is  4.7  m 
at  290  s  (or  about  1800  m),  the  maximum  depth  of 
the  probe.  We  tend  to  prefer  the  simpler  quadratic. 

4.  Concluding  remarks 

We  have  analyzed  34  CTD-T-5  XBT  pairs  of  profiles 
from  the  Sargasso  Sea  to  assess  the  temperature  and 
fall-rate  accuracies  of  the  probes.  We  found  tempera¬ 
ture  to  be  biased  about  0.07°C  (XBT  warmer),  with 
some  suggestion  that  the  offset  might  be  acquisition- 
system  specific.  A  similar  offset  was  found  by  Kennedy 
et  al.  ( 1989)  and  Wright  and  Szabados  ( 1989).  With 
one  of  the  three  datasets,  the  offset  appeared  to  be  tem- 


T-vble  4.  Summary  uf  the  be.st  two  fall-rale  equations  found  in 
this  study.  The  cubic  model  is  slightly  better  in  a  statistical  sense. 
The  accuracy  is  taken  as  the  standard  deviation  of  the  residuals 
resulting  from  the  regression  calculations  (Figs.  .Sa  and  Sbl. 


Cuelfk'ient 

Estimate 

95'i  contidcncc  interval 

Model: 

r  ht  -t  tty  accuracy 

:  standard  deviation  of  8.4  m. 

h 

6.70.5 

6.685.  6.726 

( 

1.619  X  10  ’ 

1.721  -  10  \  1.516  10  ' 

.Model;  r  - 

a  *  hi  ^  if  V  dt\  accuracy:  standard  deviation  ot'K.4  m. 

a 

'  1  .XO.t 

.5,710.  O.IO.tS 

h 

6.795 

6.716.  6,873 

2.475  /  10  ’ 

3.231  X  10  1,719  -  10  ’ 

d 

2,148  X  10  '• 

0,210  X  10  7  4,086  X  10  ■ 

perature  (or  possibly  pressure)  dependent,  but  we  do 
not  consider  our  results  conclusive.  Once  the  offset  was 
removed,  the  temperatures  had  a  standard  deviation 
of  about  0.08°C.  which  lies  within  the  manufacturer's 
claimed  accuracy. 

The  manufacturer's  fall-rate  equation. 

r  =  6.828/  -  0.00182/-. 

was  found  to  he  accurate  to  within  ±2'J  or  ±5  m  over 
all  portions  of  the  profile,  as  claimed  (  Fig,  4).  but  two 
significantly  better  fall-rate  equations  were  calculated. 
The  cubic  equation 

r  =  -1.803  +  6.795/  -  0.002475/'  +  2.148  x  10  '7  ' 

was  slightly  better  in  a  statistical  sense  than  the  qua¬ 
dratic  forced  through  zero. 

c  =  6.705/  -  0.0()1619/x 

although  the  magnitude  of  the  improvement  was  small 
and  may  not  justify  the  additional  comple.xity.  The 
errors  in  these  equations  are  nearly  depth  independent 
except  quite  near  the  surface,  with  a  standard  deviation 
of  8.4  m  for  both.  The  significant  improvement  re¬ 
sulting  from  the  use  of  these  corrections  is  illustrated 
in  Figs.  6a  and  6b. 

Several  recommendations  and  observations  come 
out  of  this  study.  The  question  of  whether  or  not  there 
is  a  small  offset  in  the  T-5  temperatures  should  be  re¬ 
solved.  and  if  so,  if  it  is  temperature  or  pre.ssure  de¬ 
pendent.  If  there  is  an  offset,  ts  it  in  the  thermistors 
themselves  or  in  the  elcctnvnics  of  the  acquisition  sys¬ 
tems?  If  there  is  a  bias  originating  in  the  thermistors 
themselves,  is  it  lot  dependent? 
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Fig.  6.  Mean  lemperature  ditTerence  versus  depth  for  the  CTD- 
XBT  profile  pairs  in  this  study.  Panel  (a)  was  computed  using  un- 
corrccled  temperatures  and  the  manufacturer’s  fall-rate  equation. 
Panel  ( b )  was  computed  using  the  temperature  corrections  and  cubic 
fall-rale  equation  from  this  study. 


The  question  of  whether  or  not  one  fall-rate  equation 
is  suitable  for  all  parts  of  the  World  Ocean  has  not 
been  answered.  Seaver  and  Kuleshov  ( 1 982 )  point  out 
the  kinematic  viscosity  of  water  changes  by  42%  be¬ 
tween  10°  and  25°C,  This  changes  the  value  of  the 
bulk  drag  coefficient,  which,  as  the  a.nalyses  of  Green 
(1984)  and  Hallock  and  T eague  ( 1 992 )  show,  is  a  ma¬ 
jor  parameter  controlling  the  fall  rate  of  the  probe.  Will 
the  different  temperature  structure  in  various  parts  of 
the  World  Ocean  lead  to  significantly  different  fall-rate 
equations?  Several  years  ago  the  authors  examined  fall- 
rate  corrections  for  T-5  XBTs  in  the  northeast  Pacific 
and  found  a  different  fall-rate  equation  for  Sippican 
T-5  XBTs.  It  was  derived  from  a  dataset  of  poorer 
quality  than  the  one  in  this  study,  so  we  are  uncom¬ 
fortable  asserting  it  is  necessarily  more  appropriate  for 
the  Pacific,  and  our  present  equations  are  appropriate 
for  the  Atlantic  (or  Sargasso  Sea,  anyway).  However, 
we  feel  the  question  is  still  open  as  to  whether  a  geo¬ 
graphically  universal  fall-rate  equation  can  be  derived 
for  each  type  of  expendable  probe. 

Both  questions  of  temperature  offsets  and  universal 
fall-rate  equations  can  probably  be  answered  only  by 
dedicated  efforts  with  sufficiently  large  numbers  of  ex¬ 
pendables,  carefully  calibrated  CTD  systems,  and  vir¬ 
tual  simultaneity  between  each  XBT  drop  and  CTD 
cast.  Most  studies — including  this  one — have  used  data 
that  were  collected  during  experiments  that  did  not 
have  XBT  calibration  as  their  primary  purpose.  Be¬ 
cause  of  this,  compromises  had  to  be  made  that  nec¬ 
essarily  qualified  the  amount  and  quality  of  the  data 


obtained.  With  the  enormous  impact  XB  I  s  have  upon 
international  databases,  it  seems  appropriate  that  na¬ 
tional  funding  agencies  consider  sponsoring  delinilixe 
studies  of  these  instruments. 

Our  final  observation,  however,  is  that  the  entry  into 
the  oceanographic  market  of  a  number  of  suppliers  of 
expendable  probes  and  of  data  acquisition  systems  may 
lower  the  price  of  the  equipment,  but  it  will  complicate 
the  problem  of  devising  accurate  temperature  and  fall- 
rate  conversion  equations  to  maximize  XBT  accuraev  . 
Each  manufacturer’s  system  and  each  manufacturer’s 
probe  has  the  potential  of  requiring  different  conversion 
equations.  Furthermore,  as  the  manufacturers  modify 
their  systems  and  probes  so  as  to  improve  their  com¬ 
petitive  position  in  the  marketplace,  new  conversion 
equations  may  be  required.  The  authors  have  first-hand 
experience  with  these  problems  with  air-deployed 
XBTs  (AXBTs);  there  is  every  reason  to  fear  the 
oceanographic  community  will  have  to  face  this  prob¬ 
lem  with  the  ship-deployed  instruments  as  well.  Indi¬ 
vidual  research  projects  may  circumvent  these  prob¬ 
lems  by  doing  an  XBT  calibration  study  each  time 
significant  numbers  of  XBTs  are  deployed  or  when 
high  accuracies  are  required,  but  this  will  not  solve  the 
problem  of  how  to  improve  the  accuracy  of  the  va.st 
amounts  of  XBT  data  being  sent  to  international  ar¬ 
chives. 
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